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Low-frequency mangnetosonic pulses in a two-ion-species plasma arc studied theoretically and 
by simulation with a one-dimensional electromagnetic simulation code based on a three-fluid 
model, with particular attention to the dynamics of minority heavy ions. It is found that 
heavy ions can gain some energy from the pulses. Because of this energy transfer, the pulses 
are damped even if the plasma is collisionless and pulse propagation is perpendicular to the 
magnetic field. 
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The presence of multiple ion species introduces many interesting effects to linear and nonlinear 
magnetosonic waves. 1-9 ) First, the magnetosonic wave in a two-ion-species plasma is split into two 
modes: the high- and low-frequency modes.& The frequency of the low-frequency mode tends to 
zero as the wave number k tends to zero and approaches the ion-ion hybrid resonance frequency as 
k — > oo.i) The high-frequency mode has a cut-off frequency on the order of ion cyclotron frequency 
and has a resonance frequency on the order of the lower hybrid frequency. (The dispersion curves 
can be found in ref .1). It was found that the nonlinear behavior of these modes are both described 
by the Korteweg-de Vries (KdV) equation, even though their dispersion curves are quite different 
in the long-wavelength region.^ The characteristic soliton widths are the electron skin depth for 
the high-frequency mode and the ion skin depth for the low-frequency mode. 

In a collisionless, single- ion-species plasma, small- amplitude magnetosonic waves propagating 
perpendicular to a magnetic field are not damped.00) (Large-amplitude waves can be damped, 
because they accelerate a fraction of the ions by the electrostatic field. 12 ~ 15 )) It was reported in 
refs. 4 and 5, however, that a nonlinear pulse of the high-frequency mode is damped in a multi- 
ion-species plasma. (Here, light ions are assumed to be the main component, as in space plasmas.) 
The damping is due to the energy transfer to heavy ions. That is, the transverse electric field in 
a high-frequency-mode pulse can accelerate heavy ions in the direction parallel to the wavefront.& 
Thus the kinetic energy of heavy ions is increased behind the pulse. This could be an important 



l 



dissipation mechanism in a collisionless multi-ion-species plasma such as the solar corona. 

In this letter, we will show that the low- frequency mode pulse can also impart some energy to 
heavy ions in the same way. Hence the pulse can be damped. An important difference is, however, 
that the heavy ions can suffer cyclotron oscillation a few times while they are in the pulse region, 
because the soliton width of the low-frequency mode is quite large, about (rn-i/m e ) 1//2 times as large 
as that of the high-frequency mode, where m; is the ion mass and m e is the electron mass. The 
heavy ions therefore alternately gain and lose energies in the pulse region. Here, we will obtain 
the net change in the kinetic energy of the heavy ions and study the wave damping due to this 
mechanism. First, we will analytically discuss the heavy- ion motion in the low- frequency mode 
pulse. It is found that heavy ions can be slightly accelerated by the pulse. Next, using a simulation 
based on a three-fluid model, we will show that the low-frequency-mode pulse imparts some energy 
to heavy ions as theoretically predicted. Therefore, the pulse is gradually damped. 

We consider magnetosonic waves propagating perpendicular to a magnetic field in a plasma 
containing two ion species; they are denoted by a and b, and the ion cyclotron frequency of ions 
b is assumed to be lower than that of a (0& < O a ). In a two-ion-species plasma, the dispersion 
relation of the low-frequency mode is given by 

u = v p0 k(l - k 2 d 2 /2) (1) 

in the long-wavelength region. Here v p q is defined as 

v p0 = v A /(l + vl/c 2 ) 1 / 2 , (2) 

where va is the Alfven speed, va = Bo/^irpo) 1 / 2 with po the average mass density po = n a Qm a + 
ribomb. The length d in eq. @) is defined as 
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where £l e is negative and uj p j is the plasma frequency of particle species j. 

When two ion species are present at considerable densities, the first term in the square brackets, 
which is proportional to 1 — f^ 1 ) 2 , is the dominant term, and d is on the order of the ion skin 
depth, c/cjpi. The dispersion of the low- frequency mode is ~ (rrii/m e ) times as large as that of the 
magnetosonic wave in a single-ion-species plasmaM In the limit of ra&o — > 0, the first and second 
terms disappear, and the electron inertia effect appearing in the third term becomes important. 
The displacement currents produce the fourth term J2j w pj/^j' which is about Q 2 e /u pe times as 
large as the third term. 

As can be expected from eq. (|l]), the nonlinear low-frequency wave can be described by the KdV 
equation, 
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where £ and r are stretched coordinates 



r = e 3 /\ 



(5) 
(6) 



with e being the smallness parameter on the order of the amplitude of, for instance, the magnetic 
field: 



B Z = B Q + eBi + e 2 B 2 + 



(7) 



Here, it was assumed that the waves propagate in the x direction and the magnetic field is in the 
z direction. The magnetic field profile of the solitary wave can be given by 

' x — Mvpot" 



B z {x,t) = B 



1 + B n sech z 
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Here, B n is the normalized wave amplitude, B n = cBi/Bq, D is the soliton width, 

D = 2(v A /v p0 )dB- 1 / 2 , 
and M is the Mach number related to the wave amplitude B n through 

M = 1 + (v p0 /v A ) 2 B Q /2. 
The longitudinal electric field E x and the transverse electric field E y have profiles 
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If rib is taken to be zero, then the above KdV theory is reduced to the one for a single-ion-species 
plasma. 

We now discuss single-particle motion of heavy ions in the solitary pulse using the equations of 
motion: 

(13) 



bx 



dt 



E x (x,t) + ^B z (x,t) 
c 



dv by 
dt 



E y (x,t)-^B z (x,t) 



(14) 



We note that the relation E y — v x B z /c = holds for both ions and electrons in nonlinear mag- 
netosonic waves in a single-ion-species plasma.lll0^ Thus E y does not accelerate ions in such a 
case. The KdV theory for the low-frequency mode also gives the relation v bx = cE y /Bo among the 
lowest-order perturbations^ The discussion in ref. 3, however, shows that this relation can easily 
break down for heavy ions. Here we calculate, more accurately, the heavy-ion motion, including 
inertial effects, and examine if the heavy ions can gain energies. 
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We integrate eqs. ( |l3| ) and (14) over time, assuming that the initial particle velocity is zero, 
v bx(fy = v by(fy = 0; an d the initial particle position is in the far upstream region, x(0) = xq~^> D. 
We disregard the perturbation of the magnetic field and approximate the variable x in E x and E y 
as x = xq. We take the time derivative of eq. ([H]) and eliminate v bx using eq. (|l3[), and obtain 



(15) 



Here F(t) is defined as 



F(t) = ^ 
v A d 
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We apply the Laplace transform to eq. (p!5[). Then, with the aid of the inversion formula, the 
velocity Vb y for t > can be found as 



Vby(t) 
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sin[f2;,(t — u)]F(u)du. 
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We can also have Vb x as 
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To obtain the velocity Vb y in the downstream region at large t, we assume that (xq — Mv p ot)/D 
— oo. Then, Vb y can be obtained as 
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where Vb m is defined as 
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If we neglect the effects of the displacement currents and the electron inertia, Vb m can be written 



as 
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Equation ( J19| ) indicates that heavy ions gyrate with the speed u& m in the downstream region. Thus, 
heavy ions have finite speed Vb m behind the pulse, even though their initial speed was assumed to 
be zero. 

Figure 1 shows Vb m as a function of the wave amplitude B n in a hydrogen-helium plasma where 
the mass and charge ratios are mb/m a = 4 and q b /q a = 2, respectively. The solid line is for the 
plasma with the density ratio rib/n a = 0.1 and the dotted line is for the plasma with nb/n a = 0.05. 
The speed Vb m increases with the wave amplitude B n and decreases with the density ratio n&/n a . 



4 




Fig. 1. Heavy-ion speed behind the pulse, U(, m , as a function of the wave amplitude B n . The solid and dotted lines 
are for the plasmas with nb/n a =0.1 and rib/n a — 0.05, respectively. 



If we substitute the soliton solution for the high-frequency modern' in eqs. (UD and (jl|), the 
heavy-ion speed behind the high-frequency-mode pulse is obtained as 



Vbm = 4nv h a'-£ „ pc p ° „ 1 - — cosech — . (22) 

Ha + ^pfe) 2 V n a J \v h u pe Bi /2 (l + B n /2a')J 

The speed is slightly greater than v a and a' is an order-unity quantity; for their precise expres- 
sions, see ref. 3. In the limit that the transit time is much shorter than the ion cyclotron period, 
eq. (22) is reduced to the one derived in the previous paper, eq. (A. 14) in ref. 3. 

Equation (jl^) gives v by {t) in the far downstream region. To observe the time variation of Vb y 
in and outside the pulse region, we numerically integrated eqs. (|l^) and (|14|), substituting the 
low-frequency-mode pulse given by eqs. (p|)-(|T2|). The amplitude of the pulse was taken to be 
B n = 0.03. Figure 2 shows the time variation of Vb y thus obtained. When a heavy ion enters the 
pulse at time VL^t = 15, Vb y starts to increase. Because the Lorentz force converts Vb y to v bx , Vb y 
decreases after reaching the maximum value. At about Vt^t = 25, the particle goes out of the pulse 
and starts to gyrate with the speed Vb m = 6.0 x 10~ 3 ^a- 

Next, we discuss the energy change of a heavy-ion particle, 

AE b = q b [ E x v x dt + q b [ E y v y dt, (23) 
Jo Jo 

using eqs. (|l~7|) and (|l~8|). The energy gain from the low- frequency- mode pulse can be written as 
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Fig. 2. Time variation of Vb v obtained by numerical integration of the equation of motion. 



Qb J Q E x v x dt = —v bm 2^ -^^b 



EyVydt 
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(24) 



(25) 



The energy gain from the longitudinal electric field E x is negative, and that from the transverse 
electric field E y is positive. Hence, the heavy ions obtain energies from the transverse electric field. 
The net change in the energy is 

AE b = m b v 2 bm /2. (26) 

It is thus expected that the pulse of the low-frequency mode is damped due to this energy transfer. 

Let us now study the propagation of the low-frequency mode using a one-dimensional, fully 
electromagnetic code based on the three-fluid model: 



+ V • (njvj) = 0, 



ra,' 



at + {vj 



qjE + ^-vj x B, 



ldB 

c~dt 



-V x E, 
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(27) 
(28) 
(29) 
(30) 



The x component of eq. (|30| ) gives the longitudinal electric field E x . The periodic boundary 
conditions are assumed. As the initial wave profiles, we use the solitary wave solutions obtained 
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from the KdV equation for the low-frequency mode, eq. (H), and observe their evolution. First, we 
confirmed that solitary pulses are not damped in a single- ion-species plasma (n& = 0). We then 
simulated a hydrogen-helium plasma and chose the hydrogen-to-electron mass ratio as m a /m e = 
1000. The density ratio was nb/n a = 0.1, and the magnetic field strength was |r2 e |/w pc = 0.5, so 
that c/va = 68.3. 
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Fig. 3. Profiles of the magnetic field, light-ion velocity v ay , and heavy-ion velocity Vb y at Q^t = 25. The initial 
amplitude is B n — 0.03. 



In Fig. 3 we show profiles of the magnetic field, light-ion (H) velocity v ay , and heavy- ion (He) 
velocity v^ y at time Vt^t = 25. The magnetic field and velocities are normalized by Bq and i>a, 
respectively. At time t = 0, the amplitude was B n (0) = 0.03, and the center of the pulse was 
at x/(c/u pe ) ~ 1100. Figure 3 shows that the heavy ions have a finite speed behind the pulse, 
as predicted by eq. (|l~9|) and Fig. 2. This heavy-ion motion produces perturbations of other 
components behind the pulse, 1000^,x/(c/u;pe)~2300, which can be seen in the profiles of v ay and 
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B z ; even in the profile of B z , we notice small-amplitude perturbations behind the main pulse. 

The observed heavy-ion speed V{, m , 2.2 x 10 -3 t>A, is smaller than the theoretical value given by 
eq. (p0|), Vb m = 7.7 x 10~ 3 wa- In deriving eq. fl2C|), we have not included the effect of the generation 
of hydromagnetic perturbations behind the pulse. Therefore, eq. ( pO| ) is interpreted to give the 
upper limit of the velocity of accelerated heavy ions. 

Figure 4 shows the time variation of the wave energy E w (t) of the main pulse. Even though we 
observe small-amplitude fluctuation, the main pulse certainly loses energy gradually. 
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Fig. 4. Time variation of wave energy E m (t) of the main pulse. 



In summary, the propagation of the low-frequency magnetosonic pulses in a two-ion-species 
plasma was studied theoretically and by simulation with a one-dimensional electromagnetic code 
based on the three- fluid model. First, the heavy-ion motion in the nonlinear pulse of the low- 
frequency mode was theoretically discussed. It was found that heavy ions gain some energy from 
the transverse electric field formed in the pulse. Next, using the simulation, we showed that the 
solitary pulse of the low-frequency mode is gradually damped. 
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